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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC).  The  results  were  obtained  by  ARO, 
Inc.,  AEDC  Group  (a  Sverdrup  Corporation  Company)  under  ARO  Project  Number 
V31K-37.  The  Air  Force  project  manager  was  A.  E.  Dietz.  The  manuscript  was  submitted 
for  publication  on  October  9,  1980. 

Mr.  B.  P.  Curry  is  presently  employed  by  Calspan  Field  Services,  Inc.,  AEDC  Division. 
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1.0  INTRODUCTION 

A  recent  series  of  experiments  was  performed  at  Arnold  Engineering  Development 
Center  (AEDC)  to  investigate  the  feasibility  of  improving  the  quality  of  System  Generated 
Electromagnetic  Pulse  (SGEMP)  simulation  in  AEDC  space  chambers.  These  experiments, 
which  were  carried  out  in  a  12- foot  sphere  in  the  von  Karman  Gas  Dynamics  Facility 
(VKF)  at  AEDC,  have  demonstrated  that  free-space  boundary  conditions  can  be  simulated 
in  a  test  chamber  by  employing  suitably  designed  damping  grids  to  attenuate  reflections 
from  the  chamber  surface.  A  centrally  located,  pulsed  biconical  antenna  was  used  in  the 
experiments  as  both  the  source  and  the  detector  of  transient  electromagnetic  waves; 
spherically  symmetric  dampers  were  located  at  various  radial  positions  between  the  source 
and  the  sphere  surface,  as  shown  in  Fig.  1 .  The  results  of  these  experiments  and  the  damping 
theory  presented  in  this  report  will  provide  a  basis  for  SGEMP  facility  design. 

The  purposes  of  the  theoretical  treatment  developed  in  this  report  are  to  establish  a  basis 
for  comparing  existing  experimental  data  with  theoretical  predictions,  to  validate  the 
experiments  and  theory,  and  to  help  determine  optimum  damping  configurations. 

Descriptions  are  presented  of  the  EMP  propagation  formalism  and  the  methods  of 
solution  which  will  be  employed  in  a  computer  program  to  be  based  on  this  work.  The  report 
does  not  include  a  description  of  the  code. 

r 

2.0  GENERAL  PROPERTIES  OF  ELECTROMAGNETIC  FIELDS 

2.1  WAVE  PROPAGATION  IN  A  DISSIPATIVE  MEDIUM 

CONTAINING  SOURCES 

To  discuss  electromagnetic  waves  which  propagate  from  a  central  source  region,  it  is 
useful  to  first  state  Maxwell’s  equations  (in  Gaussian  units)  for  fields  arising  from  specified 
sources  embedded  in  a  medium  that  is  isotropic  everywhere  except  at  specified  interfaces  and 
boundaries  (Ref.  1): 

V  *  D  =  4 rrp  ,  Vx  II  =  4ff  —  +  -  — 

c  c  e  <?t 

V  •  B  =  0  ,  V  X  E  =  -  1^5  (1) 

c  di 

It  is  assumed  that  field  quantities  and  induced  quantities  are  related  to  each  other  linearly  as 

B  =  /i  H  ,  D  =  r  E  |  and  J  =  o-  E  +  J0  (2) 
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where  ft,  e,  and  a  are,  respectively,  the  magnetic  permeability,  dielectric  permittivity,  and 
electrical  conductivity  of  the  medium.  It  is  assumed  that  the  source  current  density,  JQ,  and 
charge  density,  £>0,  exist  only  within  a  specified  region.  In  addition,  these  quantities  must 
satisfy  the  continuity  equation 

3|>  (r'.l) 

- S- - +  V'-J0<r',t)  =  0  (3) 


whereT'  is  a  vector  from  the  origin  to  an  infinitestimal  element  of  the  source  distribution. 
Using  standard  vector  manipulations  (see,  for  example.  Refs.  1  and  2),  one  can  derive  the 
time-dependent  electric  and  magnetic  field  wave  equations  from  Eqs.  (1)  and  (2). 


V  2e-^ 


pt  / d2  It,  4rr<7(9E  \  _  4 rr p  ^ o  +  ^  VPe 
+  ~  ~  o2  +  n~ 


V 


4  1TQ 
f 


(4) 


Solution  of  these  equations  is  greatly  facilitated  by  transformation  into  the  frequency 
domain,  by  using  the  following  convention  for  Fourier  transform  pairs: 


Q  (r  ,  cj) 


0(i\  t>  eiwl 


dt 


Q  (r,  t)  =  —  f  Q  (r,  <u) e  ,£Jt  dfc>  (5) 

277  J—oa 

In  Eq.  (5),  the  functions  Q  (rf  to)  and  Q  (rt  t)  may  be  either  vectors  or  scalars,  and  the  use  of 
the  tilde  over  a  quantity  denotes  its  evaluation  in  the  frequency  domain.  The  convention 
used  here  is  that  of  most  modern  physics  texts.  Electrical  engineering  texts  and  some  optics 
texts  use  the  convention  of  replacing  “i”  with  "-i”  (sometimes.stated  as  “-j”)  in  Eq.  (5)  and 
succeeding  equations. 


The  frequency  domain  wave  equations  for  the  electric  and  magnetic  fields  are  obtained 
from  Eqs.  (3)  through  (5)  and  are  stated  as 


V2  E  +  k2  E 


L  c2  0  toe  J 


?2H  +  k2  If  =  V  xt 


(6) 
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Gradient  divergence  of  the  current  in  the  electric  Field  equation  results  from  the  occurrence 
of  charge  fluctuations  in  the  source  region;  the  gradient  divergence  causes  an  essential 
singularity  when  Eq.  (6)  is  used  to  obtain  the  fields  within  the  source  region.  We  shall  see, 
later,  that  the  singularity  can  be  removed  by  use  of  a  gauge  transformation. 

The  solutions  of  Eq.  (6)  will  be  stated  in  terms  of  Green’s  functions,  which  satisfy 
Helmholtz  equations  of  the  form  (Ref.  1) 

V2  G  (r.  r ')  +  k2  G  U,  r ')  =  -S{ r  -  r ')  (?) 

whereT  is  a  vector  from  the  origin  to  the  field  observation  point,  and"r'  is  a  vector  from  the 
origin  to  a  source  point.  These  Green’s  functions  have  the  form 

C(;,r')-  oxp(iklr  —  r'l>  (8) 

4w  | r  —  r  | 


2.2  PROPAGATION  CONSTANT,  DIELECTRIC  CONSTANT, 
AND  WAVE  IMPEDANCE 


The  propagation  wave  number,  k,  which  appears  in  Eqs.  (6)  through  (8),  is  defined  by 


4  rr  a  i 

ti) 


(9) 


where  e  is  the  complex  dielectric  constant  characteristic  of  a  dissipative  medium.  The  real 
and  imaginary  parts  of  the  propagation  constant  are  stated  (Ref.  2)  as 


GO) 


In  all  materials  except  good  conductors,  the  propagation  constant  is,  approximately,  k  = 
k0  [1  +  (2ir(Ti)/(eaj)],  where  k0  =  cj/c.  By  contrast,  in  good  conductors,  the  propagation 
constant  reduces  to  k  =  (1  +  i)/Sc,  where  6C  is  the  skin  depth  of  the  conductor,  defined  as 


S  = 


■^2  7?  pai 


(ii) 
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It  is  noteworthy  that  tan  - 1  [(4xir)/(e<«>)]  is  the  angle  by  which  displacement  current  lags 
behind  the„true  conduction  current  as  a  wave  propagates  through  a  dissipative  medium.  A 
related  concept  is  the  impedance  presented  to  the  wave  by  the  medium: 


where  Za  =  {4x/c)Vfi7e‘.  The  impedance  of  a  good  conductor  is  Z  =  (1  -  i)/(6tr).  It  is  well 
known  that  the  real  part  of  the  impedance  is  purely  resistive,  and  the  sign  of  the  imaginary 
part  (reactance)  depends  upon  whether  wave  energy  is  stored  capacitively  or  inductively.  The 
biconical  antenna  used  in  the  experiments  cited  earlier  is  predominantly  capacitive,  and  the 
dampers  are  assumed  to  be  purely  resistive.  The  impedance  of  free  space  is  Z*p  =  4x/c  = 
4.192  x  10“ 10  sec/cm  in  Gaussian  units.  [Z^p  is  equal  to  120  x  ohms  in  international  units 
(SI).]  In  all  damping  calculations  the  ratio  of  the  damper  impedance  to  the  impedance  of  the 
intervening  medium  occurs;  hence,  the  choice  of  units  is  immaterial  to  the  formulation  of 
the  multipole  expansion  boundary  value  problem. 

2.3  ANTENNA  IMPEDANCE 


The  concept  of  antenna  input  impedance  is  considerably  different  from  the  wave 
impedance  discussed  above.  It  incorporates  the  resistive  effect  related  to  the  loss  of  energy 
by  radiation.  If  the  antenna  were  sufficiently  long,  the  effect  of  the  caps  could  be  ignored  in 
computing  antenna  impedance  because  the  current  would  be  damped  (by  radiation 
resistance)  enough  to  suppress  reflection  of  the  current  pulse  from  the  antenna  caps.  Since 
this  is  not  the  case  at  frequencies  of  interest  in  the  present  experiments,  the  antenna 
impedance  must  be  computed  by  dividing  the  voltage  difference  between  the  antenna  feed 
terminals  (calculated  along  an  arc  of  constant  radius)  by  the  total  current  flowing  through  a 
circular  area  whose  radius  is  equal  to  the  antenna  terminal  radius.  Equivalently,  one  may 
take  twice  the  voltage  difference  between  the  top  antenna  and  the  ground  plane  which 
bisects  the  biconical  array  between  the  feed  terminals  (see  Fig.  2).  The  required  current  can 
be  obtained  from  the  frequency-domain  form  of  Ampere’s  law. 


C 


(13) 


The  total  current  flow  is  obtained  by  integrating  Eq.  (13),  using  Stokes’  theorem 


J .  S  d  X 


1st 


f  II  •  df 


(14) 
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where  z  is  along  the  axis  of  symmetry  and,  thus,  perpendicular  to  the  ground  plane,  L  is  the 
area  of  a  circle  whose  radius  equals  the  antenna  feed  radius,  and  the  line  integral  is  evaluated 
around  that  circle.  Since  the  current  fluxes  in  the  two  halves  of  the  antenna  are  out  of  phase, 
the  electric  field  term  in  Eq.  (14)  should  vanish  by  symmetry,  leaving 

T  =  —  •  d?  =  I  -  r  sin  9  Hi  (r,  0,  <£  =  0)  (15) 

in  I  2  c 

r  <5  i 


£ 


Figure  2.  Coordinate  system  for  the  electrical  description  of  a 
biconical  antenna. 
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where  5  sin  0C  is  the  feed  radius  of  the  antenna,  and  0C  is  the  antenna  half-cone  angle.  The 
voltage  difference  between  the  two  halves  is  just 


V 


=  2 


w/  2  -2T 

f  E 

»  /I 


df 


0=0 


2 


(r,  0,  0  =  0)  d0 


Thus,  the  antenna  impedance  is 


Z 


A 


=  V/I 


4 

c 


1 

jr  Eg  <5,  e,  <f>  = 

flc _ 


0)  d0. 


S,n  H<£  lS-  ^c’  ^  =  0) 


(16) 


(17) 


In  these  equations,  the  Fields  are  evaluated  at  £  =  0  (for  convenience)  since  there  is 
axisymmetry.  Through  the  multipole  expansion  which  will  be  developed  later,  Eq.  (17) 
contains  both  the  effect  of  the  source  characteristics  and  the  damping  and  reflecting 
properties  of  the  cavity.  In  fact,  the  antenna  impedance  for  free-space  boundary  conditions 
can  be  obtained  simply  by  setting  the  multipole  coefficients  that  represent  standing  waves  in 
the  region  between  the  source  and  the  first  damper  equal  to  zero  for  Eq.  (17). 


2.4  FREQUENCY  RESPONSE  OF  THE  ANTENNA-CAVITY  SYSTEM 


After  the  Fields  have  been  obtained  and  Eq.  (17)  has  been  evaluated,  the  frequency 
response  of  the  entire  system  (antenna  plus  chamber  with  damping  grids)  is  obtained 
(assuming  that  there  are  no  losses  in  the  transmission  line  and  that  it  is  perfectly  matched  to 
the  antenna)  from  transmission  line  theory  (Ref.  3). 

TQUtW>  =  ZaU-Zc 

T  u  ZA<^  +  zc  (18) 

in 

where  Zc  =  (Zsp/7r)  fti  cot  (0/2)  is  the  (purely  resistive)  characteristic  impedance  of  the 
antenna.  The  procedure  for  comparing  theory  and  experiments  is,  thus,  the  following: 

1.  Obtain  the  Fourier  transform  of  the  input  current  pulse  (prior  to  the  arrival  of 
the  first  reflection  from  the  cavity). 

2.  Calculate  the  electric  and  magnetic  fields  from  the  multipole  expansion  which 
will  be  derived  later.  (This  step  is  optional  and  may  be  omitted  if  only  the  sensed 
pulse  waveform  is  desired.) 

3.  Calculate  the  system  frequency  response  from  Eqs.  (17)  and  (18). 
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4.  Calculate  the  inverse  Fourier  transform  Iou1  (t)  of  the  current  measured  at  the 
antenna  terminals.  This  result  can  then  be  compared  with  the  experimental  pulse 
waveform.  In  performing  Step  (4),  one  must  account  for  contributions  to  loul  (t) 
from  the  negative  frequency  domain  of  l0Ui  (o>)  as  well  as  from  the  domain  of 
positive  frequencies. 


3.0  SOLUTION  OF  THE  ELECTROMAGNETIC  FIELD  WAVE  EQUATIONS 
BY  USE  OF  HERTZ-DEBYE  POTENTIALS 

i 

The  use  of  Hertz-Debye  potentials  yields  an  automatic  separation  of  the  solutions  of  Eq. 
(6)  into  TM  and  TE  waves.  Another  advantage  of  using  these  potentials  is  their  value  as  an 
intermediary  in  the  determination  of  all  electric  and  magnetic  field  components  from  the 
radial  inhomogeneous  solutions  of  Eq.  (6).  Formal  mathematical  properties  of  these 
potentials  are  discussed  extensively  in  Ref.  4.  These  potentials  are  introduced  here  in  a 
fashion  analogous  to  that  of  Ref.  2  (pp.  28  —  30);  however,  the  treatment  presented  here  is 
generalized  to  include  dissipative  effects. 


3.1  SOLUTION  OF  THE  HOMOGENEOUS  WAVE  EQUATIONS 

The  Hertz  vectors  (ft  t)  and  Pm  (ft  t)  are  introduced  as  particular  solutions  to  the 
source  free  form  of  Eq.  (1); 

,  - 

II.  ^  -  V  x - 

1  <■  at 


r-  t  ;  05  pm 

K,  =  -c  Vk  — — 
*•  <9i 


(19) 


Additional  particular  solutions  of  Eq.  (1)  are  obtained  by  substituting  Eq.  (19)  into  Eq.  (I), 
yielding,  for  the  electric  field: 


(20) 


and 


_ _  (3  2  p 

r 7  X _ Jt. 

1  “  C2  V 
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The  analogous  equations  for  H2  are 


*"a  -  (-  *PM\ 

dt  "  v  *  V  *  di  ) 


(21) 


and 

vTlT  =  -  B. (L  +  V  X 

2  c2Ut  <  /  di 

Wave  equations  analogous  to  Eq.  (4)  are  obtained  by  use  of  the  following  procedure: 


1.  The  terms  -  v(  V*  3PE/3,)  and  -v(V«  3PM/dt)are  added  to  the  first  Ej  equation 
and  the  first  H2  equation,  respectively.  Since  the  curl  of  the  gradient  vanishes, 
this  addition  does  not  alter  the  results  of  the  next  operation  and  is,  in  fact, 
equivalent  to  a  gauge  transformation. 

2.  The  second  operation  consists  of  taking  the  curl  of  the  first  E|  equation  and  the 
first  H2  equation,  as  well  as  the  time  derivative  of  the  second  H2  equation.  Ej, 
E2,  Hi,  and  H2  are  eliminated  from  the  resulting  four  equations,  to  obtain  the 
homogeneous  wave  equations  for  the  Hertz  vectors  PE  and  Pm  in  the  time 
domain 


The  electric  and  magnetic  fields  are  constructed  from  E  =  E]  +  E2  and  H  =  Hi  +  H2, 
where  E]  is  obtained  by  inverting  the  first  part  of  Eq.  (20)  and  E2  is  given  by  the  second  part 
of  Eq.  (19).  Similarly,  Ht  is  given  by  the  first  part  of  Eq.  (19),  and  H2  is  obtained  by 
inverting  the  first  part  of  Eq.  (21).  This  process  is  greatly  facilitated  by  Fourier 
transformation.  The  electric  and  magnetic  fields  in  the  frequency  domain  are,  thus,  derived 
from  the  corresponding  Hertz  vectors  as  stated  below: 


and 


E 


1  - 


4  77  a  i 


(23) 
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The  utility  of  the  Hertz-Debye  formalism  is  made  immediately  apparent  upon  examination 
of  these  equations.  If  the  vectors  PE  and  PM  are  assumed  to  be  radial  in  form,  then  the  curl 
of  each  vector  can  have  no  radial  component.  Hence,  the  electric  field  derived  from  PM  and 
the  magnetic  field  derived  from  PE  are  entirely  transverse.  Because  of  the  terms  involving  the 
curl  of  the  Hertz  vectors,  the  electric  field  derived  from  PE  and  the  magnetic  field  derived 
from  Pm  have  non  vanishing  radial  components.  For  this  reason,  the  fields  derived  from  the 
electric  Hertz-Debye  potentials  are  called  transverse  magnetic  (TM)  waves.  Similarly,  the 
magnetic  Hertz-Debye  potential  is  the  source  of  transverse  electric  (TE)  waves.  These 
potentials  are  defined  by 

P  E  (r ,  t)  =  r  {  r  ,  t) 

and 

PM  (r,  t)  =  r  wH  (r,  t)  (24) 

It  should  be  noted  that  the  separation  of  fields  into  TM  and  TE  waves  by  use  of  these 
potentials  is  independent  of  the  choice  of  coordinate  system. 


The  wave  equation  satisfied  by  the  Hertz-Debye  potentials  is  derived  by  use  of  the 
following  vector  identities: 

V  x  r  =  0,  V^^s)  =2VS 

and 

V  x  (r  s)  =  Vs  x  r 


«  where  s  is  any  scalar.  Equation  (22)  can  now  be  reduced  to  the  form 


where 


VSxr  =  0 


S 


c2  \0i  (  )  <Jt 


V  ”E.  m 


An  equation  of  this  type  is  satisfied  either  by  S  =  0  or  S  as  a  function  only  of  r.  The  second 
choice  is  trivial;  the  first  yields  the  scalar  wave  equations  for  the  Hertz-Debye  potentials  in 
the  time  domain 


fu_  /  ^  4»?cr\ 

c2  H'  ~T~ )  By 


-V  ffRi|i 


0 


(25) 
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for  either  the  electric  or  magnetic  Hertz-Debye  potential.  Fourier  transformation  yields  the 
corresponding  scalar  Helmholtz  equations  in  the  frequency  domain: 

V2  WE,M  +  k2  "E,  M  =  0  (26) 

3.1.1  Spherical  Coordinate  Representation  of  the  Electric  and  Magnetic 
Fields  Derived  from  the  Hertz-Debye  Potentials 

The  spherical  coordinate  decomposition  of  Eq.  (23)  is  essential  for  two  reasons:  (1)  to 
permit  the  definition  of  the  boundary  value  problem  posed  by  the  dampers  in  the  spherical 
cavity  and  (2)  to  facilitate  solution  of  the  inhomogeneous  Field  equations,  Eq.  (6).  The  latter 
is  accomplished  by  a  treatment  similar  to  that  of  Ref.  4,  using  the  properties  of  the  angular 
momentum  operator  L2,  which  is  the  angular  portion  of  the  Laplacian  operator  V2  (detailed 
in  Section  3.2).  To  facilitate  the  decomposition  of  Eq.  (23),  it  is  convenient  to  use  the 
following  vector  identities: 

Vx(vs*;)  =  2Vs  +  r*v(vs)  -  r  V2  s 
and 

V2S  =  i  hi  B 

r  dr2  r2 

where  the  scalar,  s,  is  taken  to  be  either  ^  or  ttm- 


The  electric  and  magnetic  fields  derived  from  the  Hertz-Debye  potentials  can  now  be 
stated,  quite  simply,  as 


E  = 


ilTOi 

ibl 


) 


if*  co 
c  sin  6 


+(i,ils±Y' 


d  2  r Iff.  I 


i  (X&J 

3r7TM 

s2  tWe  1 

r 

-  c 

38 

6  in  8 

3r  50  J 

(27) 


and 


L2  17 m 
r 


a  r-" 

0  if  CO 
r  c 


dfyE 

3d> 


d2  ™ 

rlrM 

drdd 


+  0  d'*E  +  I d  ryM  (28) 

r_c  38  sin  £?  dr  <90 
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Analogous  equations  could  be  stated  for  other  choices  of  geometry  by  expanding  Eq.  (22)  in 
the  appropriate  coordinate  system. 

3.1.2  Boundary  Conditions  on  the  Herlz-Debye  Potentials 

The  boundary  conditions  which  are  imposed  at  each  interface  shown  in  Fig.  3  are 
continuity  of  tangential  electric  and  magnetic  field  components,'?  x  E  and  r  x  H.  Let  us 
denote  by  r~the  approach  to  the  “j th”  interface  from  outside  that  interface  —  i.e., 
r+=  rj  +  where  J*  is  a  positive  infinitesimal.  Similarly,  let  rj_=  rj  -  f  denote  the  approach 
to  the  interface  from  inside  it.  Referring  to  Eqs.  (27)  and  (28),  we  find  that  the  following 
four  matching  conditions  must  be  satisfied  at  each  interface: 


In  a  later  section,  the  multipole  expansion  of  the  Hertz-Debye  potentials  will  be  presented, 
and  a  matrix  procedure  based  on  Eq.  (29)  will  be  used  to  relate  the  multipole  coefficients  in 
all  the  regions  shown  in  Fig.  3. 


3.2  SOLUTION  OF  THE  INHOMOGENEOUS  WAVE  EQUATIONS  BY  USE  OF 
THE  HERTZ-DEBYE  POTENTIALS  AND  SOURCE  GREEN’S  FUNCTIONS 


With  reference  to  the  vector  identity 

V2  (V  t)  ^  r  •  V2  V  +  2  V  ■  V 
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(where  V  is  either  E  or  H),  we  can  show  that  the  radial  components  of  Eq.  (6)  must 
satisfy 

V2  (E  ■  r)+  k2  (E  .  d  -  -i  “  [k|  J0  ■  r  +  (2  +7  |-)  V-  J0]  (30) 


where  ZD  =  (4ir/c)/(V^/E)  is  the  (real)  impedance  of  the  medium  (which  will  usually  be  free 
space)  outside  the  source,  and 


v2  (H  ■  r)  +  k2  (H  .  r) 


Zsp  r  ^^sin  #  ^  jo  ■  § 

sin  e  L  de  &{>  . 


(3D 


where  Zsp  =  4ir/c  is  the  impedance  of  free  space.  In  what  follows  it  is  useful  to  write  the 
Laplacian  operator  as  V  h  =  [(iAj/r2]  -I-  l/r(d2rs)/(dr2)  where  L  =  -i  r  x  v  is  the  angular 
momentum  operator.  This  operator  is  extensively  used  in  both  classical  and  quantum 
mechanical  studies,  and  its  properties  (in  the  context  of  radiation  fields)  are  discussed  in 
detail  in  Refs.  1  and  4.  The  eigenfunctions  of  L2  are  ( (f+ 1).  The  solutions  of  Eq.  (26)  are 
constructed  for  a  spherical  geometry  from 

k  r  ^  M  =  2 1  m  Ct  Yf  ( 9,  <t> )  [  C\  (k  r),  ^flk  r),  ^  (k  r)  ]  (32) 


Figure  3.  Geometry  for  electromagnetic  wave  propagation  in  an 
internally  excited  hollow  sphere  containing  a  single 
damping  region. 
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where  JCp1 )  denotes  a  set  of  multipole  expansion  coefficients  which  will  later  be  determined 
for  each  region  of  the  cavity  shown  in  Fig.  1.  The  functions  fj(kr),  fftkr),  and  ^f(kr)  are, 
respectively,  Ricatti-Hankel  functions  of  the  first  and  second  kind  and  Ricatti-Bessel 
functions,  fj(kr)  and  t^fkr)  represent,  respectively,  outwardly  and  inwardly  propagating 
spherical  waves  and  V'Kkr)  represents  standing  waves  (in  undamped  regions  of  the  cavity). 
The  asymptotic  limits  of  the  Ricatti-Hankel  and  Ricatti-Bessel  functions  are,  respectively, 

j  4Tg  ~  (k  r)  =  exp|±  i  ^kr  -  ^  (  £  +  and  J^^(kr)  =  sin^kr-^2^ 

r  H  »  r  -*  OC 

Only  the  functions  ^f(kr)  are  regular  at  the  origin;  hence,  only  these  functions  can  be  used 
within  the  source  region.  Properties  of  these  functions  are  discussed  in  Ref.  5. 


The  solutions  of  Eqs.  (30)  and  (31)  are  obtained  by  use  of  the  Green’s  function,  which 
satisfies  Eq.  (7).  These  solutions  are  then  related  to  the  corresponding  Hertz-Debye 
potentials  through  the  radial  components  of  Eqs.  (27)  and  (28).  The  results  are 


L2(Sk) 


4  77  a  i 
€0) 


g  (?,  ;od3;' 

(33) 


where  Z  is  the  wave  impedance  defined  in  Eq.  (12),  and 


L  2(n)  =  H-r 

M 


z 

Sp 


36' 


(sin  0'J 


o 


h- 


<?(jo  •  e')~ 

dtp ' 


±2  d3r' 

sin  6 ' 


(34) 


where  V'  is  the  smallest  volume  completely  enclosing  the  source  distribution  and  is  not 
necessarily  spherical,  and  the  Green’s  function  is  defined  in  Eq.  (8). 

The  potentials  in  Eqs.  (33)  and  (34)  have  the  form 

k"E.V,  -Ef,  „(■[".>>, ”)  ^  (9.*)  C/U  r>  <35) 

where  jap*}  and  (bf  j  denote  the  sets  of  source  multipole  coefficients  for  TM  waves  and  TE 
waves,  respectively. 
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The  Green’s  functions  in  Eqs.  (33)  and  (34)  are  expanded  (Ref.  1)*  in  terms  of  the 
Ricatti-Bessel  functions  and  spherical  harmonics  as 


Ep  m  (kr)  Y 


^  ^  -  *1 


(*',*')  Y™ 


{8.  0) 


(36) 


With  use  of  this  expansion  and  the  orthogonality  properties  of  the  spherical  harmonics,  the 
source  multipole  coefficients  may  be  obtained  for  TM  and  TE  waves  by  premultiplying  Eqs. 
(33)  and  (34)  by  Y,^*(0,0)  and  integrating  over  a  normalizing  sphere. 


The  result  for  the  TM  wave  source  multipole  expansion  coefficients  is 


a 


m 

t 


z 

t(e+  i) 


Q 


ro-r'+(24.r' 


(kr') 

r' 


Y™V,s6'>  d3"?' 

(37) 


where*  denotes  complex  conjugate.  Similarly,  theTE  wave  source  multipole  coefficients  are 


i  z 

9 

i  m  _  _ p 

t  ~  e<e+  i) 


fff[w 


Y(m*<  (?'0')d3i 


(38) 


Equations  (37)  and  (38)  are  valid  for  a  general  source  distribution.  For  an  axisymmetric 
source,  however,  TE  waves  vanish  and  the  TM  wave  coefficients  reduce,  outside  a  spherical 
region  completely  enclosing  the  sources,  to 


Vit(  2 £  +  l) 

f(e+ 1) 


-  kd 

z  — 


//[■■ 


(kr') 

kz  1 


_b_ 

dr' 


(,'2v'i«)] 


'  Pj,  (cos  8')  sin  66' dr' 


(39) 


•The  Green’s  function  expansion  stated  in  Ref.  1  involves  the 'spherical  Bessel  functions  Jf(kr)  and  spherical 
Hankel  functions  h^(kr).  The  Ricatti  functions  are  simply  the  products  of  these  functions  with  their  arguments. 


i 
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Obviously,  one  must  know  the  source  current  distribution  to  use  these  equations.  Equations 
(37)  and  (38)  are  presented  because  of  the  necessity  to  consider  arbitrary  source  currents 
when  satellite  components  are  irradiated  in  a  test  chamber.  In  the  next  section,  Eq.  (39)  will 
be  evaluated  for  an  uncapped  biconical  antenna.  In  addition,  rigorous  coupled  integro- 
differential  equations  will  be  derived  in  order  to  calculate  the  current  distribution 
characteristic  of  a  capped  biconical  antenna.  In  succeeding  sections,  the  cavity  boundary 
value  problem  will  be  solved  and  the  resulting  regional  multipole  coefficients  will  be  related 
to  the  source  multipole  coefficients  stated  in  Eqs.  (37)  and  (38). 

3.2.1  Source  Multipole  Expansion  for  an  Uncapped  Biconical 
Antenna  Radiating  into  an  Unbounded  Medium 

Although  its  use  in  analyzing  the  existing  test  chamber  results  requires  that  Eq.  (36) 
should  incorporate  the  current  distribution  calculated  (with  the  equations  presented  in  the 
next  section)  for  a  capped  bicone,  some  of  the  properties  of  the  source  multipole  coefficients 
can  be  demonstrated  by  evaluating  Eq.  (39)  for  a  specific  current  distribution:  Schelkunoff’s 
mode  solution  (Ref.  6)  for  an  uncapped  bicone.  For  this  case,  Eq.  (39)  can  be  evaluated  in 
closed  form.  Referring  to  Fig.  2,  we  consider  a  system  of  two  cones  generated  by  revolution 
of  the  line  segment  d*5  about  the  positive  and  negative  z-axes.  We  shall  require  that  the 
radial  gap,  26,  between  the  cones  be  small. 

The  form  assumed  for  the  current  distribution  is  purely  radial: 

J w ( r" ')  =  - - -  Z  *o'  }  (<U,  r')  [^(cos  0  -  cos  0C)  -  fi(cos  6  +  cos  0r)J 

2tr  r i  -  1 

This  form  is  a  generalization  of  the  current  distribution  stated  in  Ref.  1  for  a  thin,  linear 
antenna.  The  primary  current  mode  is  sinusoidal: 

r')  =  Io  (o)  sin  (d  -  r')  J 


The  reflected  mode  is 


I*2,U>.  r')  =  -if0(oj)—  cos  [k0  (d  -  r')J 
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and  its  contribution  to  Eq.  (39)  vanishes.  The  higher-order  modes,  which  arise  from  the 
reaction  of  the  radiated  fields  on  the  antenna,  are  stated  (from  Ref.  6)  as* 


i  i  (to) 

O 

2  n 


Zs  oo 

P  r  4?  +-  3  .1 

zc  g70  <f  +  1)  (2t,  1)  ^  + 


(kd)  + 


(kor') 


(41) 


The  reader  should  note  that  I0  (w)  is  not  the  current  measured  at  the  input  terminals.  It  is  a 
complex  quantity  which  is  related  to  the  input  current  and  voltage  by  the  following  pair  of 
equations: 

I"  za 

1  (to,  8)  =  L  (cdl  sin  (k  d)  -  i  - cos(kd) 

0  L  zc 


and 


V  (to,  5)  =  Io  (to)  [zA  sin  (kd)  -  i  Zc  cos  (kd)J 


(42) 


The  antenna  impedance  in  these  equations  is  also  obtained  from  Schelkunoff  s  mode 
solution: 


ZA  = 


■(3) 


(to,  r')/T0 


(to) 


(43) 


The  form  of  the  current  distribution  in  Eq.  (40)  can  be  exploited  to  cast  Eq.(39)  into  the 
following  useful  equation,  in  the  fashion  of  Ref.  1 . 


(44) 


•Because  of  the  choice  of  outgoing  boundary  conditions  in  Ref.  6,  fjf  +  |(kd)  appears  in  the  complementary 
current  modes,  and  there  is  no  distinction  between  k  and  k0. 
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where 

T>,0  .  £T»<.,r1 

1=1 

Note  that  only  odd  0  multipoles  occur  because  of  the  symmetry  of  the  bicone.  Insertion  of 
the  three  current  modes  stated  in  Eq.  (39)  into  Eq.  (46)  yields  (as  k5— 0) 
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(n  +  1)  (2n  +  1)  (Jrftkd) 


4  n  +  3 


n,  ? 


(45) 


where 


a 

Jn,  l  55  ko/  [  ^2n+ 1  (kor°  ^{kr')- ^  ^2n  +  1  <ko  r#>  4>t  (k  r'>]  dl 


The  remaining  integral  can  be  evaluated  by  using  the  properties  of  the  Ricatti-Bessel 
functions  and  Bessel  function  integrals  and  recursion  relations  stated  in  Ref.  5.  After  much 
tedious  mathematical  manipulation,  we  find 


^n,  £ 


2n 


TiTT  j  -KTT  ( ^  -  X)  *t™**  +  i*.'™ 


to  n  +  i  (k  d) 
(0-1)  frg(kd)  \  j— — 


2n  + 


TTT [inTTj-  *i (kd) <kd> *1. +  1  <k= d)] 


d) 
(46) 

Assuming  that  |(k2/k2)  -  1|  <<  1,  we  obtain  the  final  form  for  the. source  multipole 
coefficients  of  an  uncapped  biconical  antenna: 


/  20  +  1 

I  -  <-l)E 

P £  (COS  0^  ) 
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(47) 


22 


A  E  DC-TR -80-58 


where  xf(k,  kc,  d)  is  denoted  the  “complementary  antenna  function”  and  is  defined,  in 
accordance  with  Eq.  (46),  as 


Xr(k,LA,d)  = 


1  n  ■ 


7T  7. 


2  (n  +  1 )  h  y  P 
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(48) 


It  is  interesting  to  note  from  Eq.  (47)  that  the  contribution  of  the  complementary  current 
modes  to  the  source  multipoles  is  out  of  phase  with  the  contribution  of  the  primary  mode. 
Since  these  equations  are  closed-form  expressions,  their  evaluation  permits  a  check  on  the 
source  coefficients  obtained  from  the  numerical  solutions  of  the  more  general  capped  bicone 
integrodifferential  equations  derived  in  the  next  section. 


3.2.2  Derivation  of  the  Coupled  Integrodifferential  Equations  for  the 
Current  Distribution  on  a  Capped  Biconical  Antenna  Radiating 
into  an  Unbounded  Medium 


The  presence  of  caps  on  the  bicone  complicates  the  treatment  of  the  last  section,  in  this 
section,  the  equations  describing  the  radial  current  distribution  in  the  conical  sections  of  the 
antenna  and  the  expanding  and  contracting  current  circles  in  the  end  caps  are  derived  by 
treating  the  antenna  as  a  boundary  value  problem  and  using  a  gauge  transformation  to  treat 
the  singularities  which  occur  in  Eq.  (39)  when  the  observation  and  source  points  coincide 
(r1-~r').  Using  this  transformation  introduces  the  effect  of  an  accumulation  of  fluctuating 
charge  near  the  antenna  tips.  Conservation  of  current  (Kirchoff  s  law)  at  the  cone-cap 
juncture  prevents  charge  accumulation  at  the  caps. 

The  boundary  value  problem  is  posed  by  requiring  the  tangential  electric  fields  on  the 
antenna  surfaces  to  vanish,  assuming  that  the  surfaces  are  perfect  conductors.  Referring  to 
Fig.  2,  we  can  see  that  the  vanishing  of  the  electric  field  on  the  conical  parts  of  the  antenna  is 
expressed  by 

0  x  E  =  ?  E^  -  <i  Ef  =  -  <£  Er  =  0  at  6  =  8C  and  n  -  6C,  and  for  S  <  r  <  4, 

at  0  =  6C  and  t  —  6C,  and  for  6  <  r  s  d,  where  25  cos  0C  is  the  axial  separation  distance 
between  the  cones.  The  azimuthal  field,  E^,  vanishes  by  symmetry.  The  radial  field,  Er, 
contains  contributions  from  boih  the  cones  and  the  end  caps. 
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Similarly,  the  tangential  field  on  the  end  caps  is  given  by 

z  X  E  =  ^Ep-pE^  =  <£Ep 

where  e,  0,  and  0  are  cylindrical  coordinates  describing  the  antenna,  and  the  tangential  field 
is  evaluated  at  z  =  ±  h  for  o  s  g  qq.  This  field  also  contains  contributions  from  sources 
on  both  the  cones  and  end  caps. 

The  solution  of  the  first  part  of  Eq.  (6)  is  obtained  by  using  properly  selected  singular 
Green’s  functions.  Caution  is  required,  however,  because  of  the  gradient  of  the  divergence 
of  the  current  density  in  the  first  part  of  Eq.  (6).  If  this  term  appears  inside  the  source 
integral,  as  in  Eq.  (39),  there  results  a  non-integrable  singularity  (Ref.  7)  when  observation 
and  source  points  coincide;  consequently,  the  treatment  of  Section  3.2.0  cannot  be  used  to 
describe  the  fields  on  the  antenna  surface.*  We  shall  instead  transform  the  field  in  such  a 
way  that  the  divergence  of  the  current  density  appears  inside  the  integral  and  the  gradient 
appears  outside  (as  in  Ref.  7).  Such  a  transformation  is  equivalent  to  the  specification  of  the 
gauge  condition  for  the  conventional  vector  and  scalar  potentials. 

Let  the  transformation  be 

E  =  U  +  V  0  (49) 

and  choose  \p  -  (V  •  U)/k2;  then  the  appropriate  Helmholtz  equation  is 

V2II  -t-  k2U  =  -  4ffi  —  T„  (50) 

C2 


with  the  formal  solution 


u  =  L  Z  k  fff  J0  (r  ')  G  (t,  r  ')  d3  r  '  (51) 

V' 

where  Z  is  the  wave  impedance  defined  in  Eq.  (12)  and 


G(r,  f  ')  = 


i  k  I  r 


4  ff  r  ■ 


•The  author  acknowledges  Dr.  C.  M.  Butler  of  the  University  of  Mississippi  for  bringing  this  point  to  his 
attention.  Dr.  Butler  also  supplied  a  derivation  of  the  coupled  integral  differential  equations  based  on  vector  and 
scalar  potentials. 
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is  a  Green’s  function  which  has  an  integrable  singularity  at  "r  =  "r ' ,  and  T'  and  T  are, 
respectively,  vectors  from  the  origin  to  a  current  element  and  to  the  observation  point 
(r,  9,  <t>). 

It  is  convenient,  at  this  point,  to  state  the  radial  unit  vector  of  a  source  point  in  terms  of 
the  coordinates  of  the  observation  point  for  both  spherical  and  cylindrical  coordinate 
systems  (see  Fig.  2).  For  a  source  point  on  the  surface  of  either  cone,  we  find 

*  r  '  =  p  sin  8'  cos(0  -  <jj')  +  z  cos  9'  <j>  sin  {<}>  -  <}>')  sin  6 '  (52) 

and 

r  '  =  p[  sin  0 '  sin  $  cos  (cf>  -  <f>')  +  cos  6 '  cos  (?] 

-rtHsin^'coB^eosl^-i^'J-cosf^sin^l-^sin  0^sin(^-0'|  (53) 


with 


r-0c,  rr-d 


c 


Similarly,  for  a  source  point  on  the  surface  of  either  end  cap,  we  find 


p  ’  =  p  cos(i£  -  <j> ')  -  0  sin(<£  -  0') 


(54) 


and 


p  '  -  r  sin  6  cos(0  -  <p' )  +■  0  cos  8  cos(0  -$')-$  sin^S 


(55) 


with 


*'  =  ±h 


These  projections  will  be  used  to  determine  the  magnitude  of  the  vector,  r  -  r ' ,  from  source 
point  to  observation  point  for  four  specific  cases  and  the  associated  Green’s  functions. 


The  Green’s  function  for  current  elements  and  observation  point  both  on  the  conical 
surfaces  is  now  stated  as 


')  = 


i  k  R 

«• 


I 


477  K 


with 


«t 


+  r -  2rr'  [a in  8  sin  9'  cos (^6  -  <£')  +  cos  8  cos  6^\ 

IT"-*  6  .  n  -  8 
c  c 


(56) 
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The  Green’s  function  for  current  elements  on  the  end  caps  and  observation  point  on  the 
conical  surfaces  is 


Cf,  T  ') 


ik  B2 

in  Rj 


with 


R2  =  |  \/r2+z''2  +  p'2-2r(8in  dp'  cos(0  -</>')+  z  "  cos  6) 


(57) 


z  ±  h 


The  Green’s  function  for  current  elements  and  observation  point  both  on  the  end  caps  is 

it  p 

=  -1 - 

K  4?rpi 


with 


yj p2  +  2  2  +  p'2  +  z  '2  -  2 [pp'  cos  ($  -  <f>‘)  -  z  z  "] 
z'  ->  ±  h 


(58) 


Finally,  the  Green’s  function  for  current  elements  on  the  conical  surfaces  and  observation 
point  on  the  end  caps  is 


G 


2 

be 


(  r  ,  r  ') 


ikp2 

e 

P2 


with 


P  2  - 


d'+dQ,n-6c 


2r'[sin  9'p  cos  (i£  —  <f>')  +  z  cos  9'] 


(59) 


Having  constructed  the  Green’s  function  for  the  four  cases  of  interest,  we  may  now 
determine  the  current  flowing  in  the  conical  surfaces  and  in  the  end  caps.  The  conical  source 
is  a  loop  of  current  flowing  in  the  (outward  or  inward)  radial  direction  whose  magnitude  is 
given  by 

Tbc(r')  =  2?7  sin  9C  tbc  r  '  *  Jo  (r  (®®) 
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where  tbc  is  the  thickness  of  the  very  thin  surface  layer  in  which  current  flow  occurs.  In 
obtaining  Gq.  (43),  we  have  used  the  approximation  that 


l 


6  +  (  A  07 2) 


sin  0'dff'  =  sin  0 


be 


e,  -  ( A  072) 


with 


A  0'  = 


* 

r 


The  volume  element  is,  similarly,  reduced  to 

i 

d3  r  '  =  r'tbcsin  0C  {§(0'-  0C).  S[0'- ( *  -  0C)]  }  d0'd0'dr' 

Allowing  for  the  phase  difference  between  currents  flowing  in  the  two  cones,  we  obtain, 
finally, 

-s 

1*W 

J  jjc  ( r  ')  d3  r*  '  =  rTbc(r'}  |  S(0'- 0C  )-$  [0'  -  ( ir  -  0C )  ]  }  dr'd0'd£'  (61) 

By  analogous  reasoning,  we  find  the  cap  current  density  volume  element  product  to  be 

d3  r  '  =  y  Icp(p)  [s(z'~  h)-  S(  z '  +  h)]dp'dz  '  (62) 


with 

Icp(P  =  0)  =  0 

The  contribution  of  both  cone  currents  and  cap  currents  to  the  auxiliary  function  U, 
evaluated  near  a  cone,  is  found  [from  Eqs.  (51),  (61),  and  (62)]  to  be 


U 


be 


(?,r',0')-X 


1- 

be 


(f,r',  0')]  dr' 


2rr 


f  f  P%p(p')[XcJ  (^0'.  ^')-^c2p  ^r, /»'.«')]  dp'  ^  }  (63) 
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where 


,0  +(A072) 


X,  (r,  r i£'l  = 
be 


/  -i/ 

'  Gbc  [s(0'-0c).5(0'-U-0c))]  <W' 

0C-(A072) 


and 


X 


±  [K  +  (tcp/2)l 


± 1  h  -  (‘cp/2)l 


Gc2  (?,?')  [a(Z'-h)tfl(Z'+h)]dZ* 


where  tcp  is  the  thickness  of  the  current-carrying  layer  in  the  end  cap  surfaces. 
Similarly,  the  form  of  U  evaluated  near  a  cap  is 


UCP  = 


jf  £  f  A.  Crl  [xL+ ('.  *->  ~  X  L"  <?•  *'>>'  £ 

f  "f  °  ?  [xiP+<  '•  o'-  *•>  -  xip-  ( Jo'  | 


where 


(64) 


0  +  *L' 
c  2 


*Lr<r’r-*') 


y*  ^  Ggc  (m'J  [s(0'-0c),S(0'-(ir-0c))]  d0' 


2 


(64a) 


and 


x  p',0')  E 
cp  -#+  h  _  cp 


x 


±  h  + 


cp 


cip  f  ')[«Z'-  h),  fi(Z'+  h )]  d0' 


(64b) 


In  forming  the  divergences  of  these  functions,  we  note  that  the  gradient  operates,  initially, 
only  on  the  unprimed  coordinates  in  the  Green’s  functions.  Thus,  the  divergences  have  the 
-form 
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V  *  U  =  i  Z  k  I  |  [  (n')  n*  V  (r,  n',  ^')  dn' 

J0  J  a  2Tr 


=  iZk 


-2JT  ^ 

. }. 


I  (iT) 


<9n ' 


<t>' 


dn'^ 

2  17 


(65) 


where  n'  is  the  unit  vector  in  the  direction  of  current  flow  and  n'  is  the  associated 
coordinate. 


The  latter  form  can  be  integrated  by  parts,  yielding 


V*  U’= 


tr,  n',  (ft ') 


dl(n') 

5ny 


tr,  n'9 


2n 


dn 


•} 


(66) 


In  an  ordinary  boundary  value  problem,  the  Green’s  function  would  normally  be  defined  so 
as  to  assure  vanishing  of  integrated  parts  containing  either  the  current  (Dirichlet  boundary 
conditions)  or  its  normal  derivative  (Neumann  boundary  conditions)  on  the  surface 
enclosing  the  source  elements.  In  this  case,  however,  the  boundary  conditions  which  must  be 
satisfied  are  the  vanishing  of  the  tangential  fields  at  the  surfaces  and  the  KirchofFs  law  for 
continuity  of  current  at  the  cone-cap  junctures.  The  latter  is  stated  as 

+  s  0  (67) 

Since  it  can  be  shown  that  xfp  (<".  Go.  <t>')  =  xi*(r»  d,  one  can  write  the  divergence  of 
the  auxiliary  function  near  the  conical  surfaces  [incorporating  Eq.  (67)  into  Eq.  (66)]  as 


V 


U 


be 


(r) 


-iZH 


he  <S> 


/  [*bc 


(r,  5,  61  -  (r,  *, 


* 

2n 


aTcp 

dp' 


( r ,  p',  (f> ')  - 


d  <f>' 
2  n 


dp' 


+ 


-2rr 

J  [X'bc  (r,  r',  <£')“  X1  (r. 


dr' 

2n 


} 


(68) 
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Again,  through  an  analogous  reasoning  process,  we  obtain 


■  -‘Zk^bcWjT  [xh20+(r,  5,  Ci'J-xL  ir,  5.*')] 


/p°  di  (pi  r2n  r  ,  _  ,  -i  j  j  ' 

~sr~f.  [4' m*-'- *!,- 


rd  dlbcu/)  f2,rr  2+r  ,  ,,,  2-  .  ,  ,vi  d6*  . " 

f  - r~ ; -  I  Xbc  (r,  r  ,  «^>  )-  Xt  (r,  r  ,  <t>  ]  -  dr  - 

•Jg  d  r  Ja  L  J  2tt 


(69) 


These  divergences  act  as  scalar  potentials,  as  can  be  seen  from  Eq.  (49). 


The  final  expressions  for  the  tangential  fields  are  obtained  from  Eqs.  (49),  (63),  (64), 
(68),  and  (69),  using  the  following  cross  products: 

0  x  r  '  =  -  ©  [sin  6 '  sin  6  cos  (p  -  dtl  +  cos  6'  cos  d\  -  r  sin  8  ’  sin  -  <£  1 


Q  \  p  ’  -  -  p  sin  6  cos  (p  -  d> ')  -  r  sin  (p  -  <£ ') 

9  x  V 


^  <9__  _  ± 

f)r  r  sin  9  d({> 


/.  x.  r  '  =  ci  sin  6'  cos  (<£  -  d)l  +  p  sin  9‘  sin  (0  -  <f>l 


z  x  p  '  = 


p  cos  (d>  -  (f)')  +  p  sin  (p  -  cH 


z  v  V 


g  A 

P 


(70) 


which  were  obtained  from  Eqs.  (35)  through  (38).  When  Eq.  (70)  is  used  to  obtain  the 
tangential  field  components  from  Eqs.  (49),  (63),  (64),  (68),  and  (69),  the  terms  proportional 
to  sin  (& -</>')  do  not  survive  integration,  since  these  are  odd  terms,  whereas  all  other  terms 
in  the  integrands  having  the  form  of  Eq.  (70)  are  even  terms.  In  addition,  the  aximuthal 
derivatives  vanish  by  symmetry.  Thus,  the  tangential  electric  fields  are  purely  azimuthal, 
both  near  the  cones  and  near  the  caps. 
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It  is  convenient  to  introduce  three  classes  of  kernel  functions  defined  as  follows: 


X  (r,  n') 
£  (r,  n ') 
)i(rtn') 


f  6')- x~  ^ 

Jo  2n 

Jcos  6 '  [  x+  (r  •  n  d> ')  -  X~  ( r »  "  '*  <t>  "1] 

o  2ff 

f  lx+  <r.  n',  </>'}-  x~  (r,  0')1  — 

Jo  2ff 


(71) 


We  state  the  final  field  equations,  in  terms  of  these  functions,  noting  that  cylindrical 
symmetry  permits  <£ '  to  be  set  equal  to  zero  and  with  the  order  of  integration  interchanged. 

<?xEbc  =  -i  </>  Z  k  j  sin  0c  sin  6  J*  Ibc  (r')£  bc  (  r,  r')  dr'---  cos  Sc  cos  dj  Ibc  (r ')  Dll1  (  r ,  r ')  dr ' 

tBinScJ  ”Tcp  tp") KJp (f. p’)ip’- [Tbc(s)^- Kj0(?, s) 


3  r‘  d,i.  <r')  vi  .  a  t 

rs-j,  l 


Po  •‘‘"(pi  v!  ,  .ll 
Kp(’-p‘ip  J  J 


(72) 


and 


zxEbc  =  i 


i  Z  k|  ein  ec  JT  Tbc(f')f^(r.r-)dr^j'’°rcp  (p'>£  Jp  G,  ,Mp- 

[rbe(s,|KL(;.S)v|f 


d  (r'J 


•iC^1  '»"-«•]) 


(73) 


The  coupled  integrodifferential  equations  for  the  currents,  Ibc,  on  the  cones  and  Icp  on  the 
caps,  are  obtained  by  setting  Eqs.  (72)  and  (73)  equal  to  zero  in  the  limits  0  —  0C,  r  -  0C,  and 
z  —  ±h.  The  results  are  stated  as  follows: 
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Jpd  ^  m  <J 

Ihc  (r')£hc  (r.  «■')  dr'+cos2  6C  cos  Q  f  E  (r')  HI1  <  r ,  r')  dr' 
S  JS 

+  sin  dc  j  Tcp  (p')K2p(?,  p')dp'~  [Tbc(.S)  ~  Kjc(r,  5) 

•  Hi 1 5P-  >1 >‘K(f  ^ 


(74) 


and 


8in0cf  Tbc(r')SL(^r')dr^r  °TcP  (p')£ip(r,P')dp' 

“'S  Jo 

-[T*<sl*  K*<?,4) 


*  iif  ^ 

OT  \%/g  dr 


( r,  r')  di 


aTcp(p') 

dp' 


0 

(75) 


In  using  these  equations,  note  that  t  -  (ge)  +  (z  z),  in  order  to  evaluate  terms  in  3/3q. 

There  are  several  options  for  solving  these  equations.  Conventionally,  the  method  of 
moments  (Ref.  8)  can  be  used  to  reduce  the  equations  to  a  system  of  algebraic  equations, 
provided  that  the  singularities  in  the  kernels  X,  and  )Ii  can  be  isolated.  This  method 
involves  establishing  a  computational  mesh  on  both  the  cone  surfaces  and  on  the  end  caps. 
The  currents  are  expanded  in  a  series  of  functions  chosen  for  computational  convenience 
and  defined  within  each  interval  of  the  mesh.  Weighting  functions  (known  as  testing 
functions)  are  similarly  chosen,  and  the  expanded  Eqs.  (74)  and  (75)  are  premultiplied  by 
these  functions.  Integrating  the  resulting  equations  over  the  finite  elements  of  the 
computational  mesh  then  yields  algebraic  relations  for  the  current  expansion  coefficients 
which  can  be  solved  by  matrix  methods.  A  large  body  of  literature  exists  that  centers  on  this 
method  and  the  properties  of  various  expansion  and  testing  functions  (in  the  context  of 
specific  functions)  (e.g..  Refs.  9  through  1 1).  When  the  expansion  and  testing  functions  are 
the  same,  the  procedure  is  known  as  Galerkin’s  method  (e.g..  Ref.  12).  In  view  of  the 
physical  significance  of  sinusoidal  currents  demonstrated  in  Section  3.2.1,  it  is  appealing  to 
use  piecewise  sinusoidal  functions  as  either  or  both  expansion  or  testing  functions,  provided 
that  the  resulting  integrals  are  evaluable.  It  should  be  noted,  however,  that  simpler  functions 
have  been  found  useful  (Ref.  10),  although  a  finer  mesh  may  be  required.  Further  discussion 
of  the  method  of  moments  is  beyond  the  scope  of  this  report.  It  may  be  possible  to  devise  a 
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simpler  scheme  to  solve  Eqs.  {74)  and  (75).  Regardless  of  how  the  solution  is  effected,  the 
resultant  current  distribution  must  then  be  inserted  into  Eq.  (39)  to  obtain  the  source 
multipole  coefficients  for  the  biconical  antenna.  If  the  method  of  moments  is  used,  then, 
obviously,  the  expansion  functions  should  be  chosen  with  ease  of  integration  of  Eq.  (39)  in 
mind. 

3.3  FORMULATION  OF  THE  CAVITY  BOUNDARY  VALUE  PROBLEM 

In  this  section  the  multipole  expansion  for  the  Hertz-Debye  potentials  is  obtained  for  all 
regions  of  the  spherical  cavity  external  to  the  source.  For  simplicity,  the  surface  of  the  sphere 
is  assumed  to  be  much  thicker  than  its  skin  depth  at  all  frequencies  of  interest.  The  multipole 
coefficients  in  each  region  of  the  cavity  (see  Figs.  3  and  5)  are  then  obtained  by  applying  the 
boundary  conditions  stated  in  Eq.  (29)  to  the  potential  expansion  shown  in  Eq.  (32).  We 
illustrate  the  procedure  by  deriving  the  multipole  coefficients,  first,  for  an  empty  sphere  and 
comparing  our  results,  in  the  limit  of  a  Hertzian  dipole  source,  with  those  of  Ref.  13.  After 
this  discussion,  the  matrix  formulation  for  N  dampers  is  presented,  and  the  results  for  one 
and  two  dampers  are  derived  in  detail.  Throughout  this  section  the  full  expansion 
corresponding  to  an  arbitrary  source  distribution  [Eqs.  (37)  and  (38)]  will  be  retained,  rather 
than  assuming  source  axisymmetry.  The  following  notation  convention  will  be  employed: 

L.  Capital  English  letters  denote  coefficients  of  standing  waves  (Ricatti- Bessel 
functions  with  real  arguments). 

2.  Lower-case  English  letters  denote  coefficients  of  outwardly  propagating  (plus 
left  superscript)  or  inwardly  propagating  (minus  left  superscript)  undamped 
spherical  waves  (Ricatti-Hankel  functions  of  the  first  and  second  kind, 
respectively,  with  real  arguments). 

3.  Greek  letters  denote  coefficients  of  outwardly  propagating  (plus  left  superscript) 
or  inwardly  propagating  (minus  left  superscript)  damped  spherical  waves 
(Ricatti-Hankel  functions  with  complex  arguments). 

To  illustrate  this  convention,. in  reference  to  Fig.  3,  the  multipole  coefficients  for  the 
single  damper  configuration  are  related  to  the  corresponding  spatial  regions  as  follows: 

a.  Source  region  -  aj>\  bf1 

b.  Region  1  -  between  source  and  damper  -  Ap1,  Bp* 

c.  Region  II  -  damper  -  ^p1,  ^p1 
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d.  Region  III  -  between  damper  and  surface  -  Ep,  Mp  ;  _ep,  _mp 

e.  Surface  region  -  +yp,  +5p 

Note  that  the  first  of  each  pair  of  coefficients  denotes  the  TM  wave  multipole  coefficient 
and  the  second  denotes  the  TE  wave  multipole  coefficient-  We  shall  return  to  the  single 
damper  problem  later,  but  now  we  shall  consider  the  empty  cavity  problem  to  illustrate  the 
formulation  and  solution  of  the  boundary  value  problem  relating  all  other  multipole 
coefficients  to  the  source  multipoles  ap  and  bp,  whose  relationships  to  the  source  current 
distribution  have  been  stated  in  Eqs.  (37)  and  (38). 

3.3.1  Multipole  Expansion  Coefficients  for  an  Empty  Cavity 

From  Eqs.  (29)  and  (30)  we  obtain  the  relations  between  inerior  and  exterior  solutions* 
for  the  TM  waves  in  a  cavity  containing  a  nondissipative  medium: 


|Ssf  [.JfJ(as).A»^(a,l]  - 


and 


ae  ^  K)  +  A?^(fl5)  Vj? 


1  +■ 


4  17  £7  i 


(76) 


where  the  subscript  “s”  denotes  surface,  as  =  k0R,  0,  =  ks/k0,  and  the  ratio  of  interior  and 
surface  impedances  is  Z/Zs  -  05  (/i/ps),  with 


—  f  < 


(4.7  o  i  \ 


1  +  i 


The  corresponding  equations  for  TE  waves  are 

and 

b™  C\\as)  +  B™  U£)  -  as) 


(77) 


♦Henceforth,  primes  denote  differentiation  of  a  function  with  respect  to  its  argument. 
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Eliminating  the  surface  multipole  coefficients,  +7jP  and  +5JP,  from  these  equations 
produces  the  cavity  multipole  coefficients  (i.e.,  the  coefficients  of  standing  waves  in  the 
cavity)  for  both  the  TM  and  TE  waves 

^  E _ _ _  {  (78) 

and 

B™  =  - 2 - - -  \  (79) 

(  f-  Z  g  <0. 


These  standing  wave  coefficients  have  an  obvious  similarity  to  the  well  known  Mie- 
scattering  coefficients — except  that  here  the  problem  is  one  of  an  empty  cavity  surrounded 
by  an  electrically  infinite,  absorbing  medium,  whereas  the  usual  Mie-scattering  problem  is 
one  of  an  absorbing  and  dispersive  sphere  embedded  in  "a  nonabsorptive  medium.  These 
coefficients  are  now  used  to  obtain  the  cavity  fields  for  comparison  with  an  earlier  analysis. 

In  an  early  study,  Butler  and  Van  Bladel  (Ref.  13)  obtained  the  fields  inside  a  dipole- 
excited  hollow  sphere  embedded  in  an  infinite  dielectric  medium.  To  compare  our  results 
with  theirs,  we  form  the  ratio  of  the  tangential  electric  and  magnetic  fields  at  the  inner 
surface  of  the  sphere: 

E0  k8C  /  j  <?airE  \ 

~  ir»w  ' dade  /  36  A  '  (80) 

since  the  azimuthal  derivatives  in  Eqs.  (27)  and  (28)  vanish  because  of  axisymmetry.  For  a 
small  dipole,  |kd|  <  <1,  only  the  first  multipole  is  important,  and,  using 


we  obtain 


where 
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This  is  the  result  obtained  by  Butler  and  Van  Bladel  (Ref.  13).  Note  that  our  choice  of  the 
sign  of  a)  [from  Eq.  (5)]  dictates  the  use  of  h^as)  for  correct  outgoing  boundary  conditions, 
whereas  Butler  and  Van  Bladel  use  the  opposite  sign  for  w  and,  hence  h^ag)  for  outgoing 
waves.  Also,  our  use  of  Gaussian  units  introduces  a  factor  C  into  Eq.  (81),  in  contrast  to 
Butler  and  Van  Bladel’s  form  of  Eq.  (83). 

3.3.2  Matrix  Formulation  for  the  multipole  Coefficients  for  the 
Multipole  Coefficients  in  a  Cavity  Containing  N  Dampers 

The  most  systematic  way  of  incorporating  the  boundary  conditions  stated  in  Eq.  (29) 
into  the  boundary  value  problem  for  a  cavity  containing  several  dampers  is  to  use  matrix 
procedures.  Let  us  define  the  source  column  vectors,  whose  dimensions  are  2(2N  +  1),  Ye 
and_Yjvii  as  vectors  whose  first  elements  are  proportional  to  af  and  bp1,  respectively,  and 
whose  remaining  elements  are  all  zero.  The  column  vectors  for  the  unknown  multipole 
coefficients  are  denoted,  accordingly,  as  Xe  and  Xm-  These  vectors  are  obtained  by 
inversion  of  the  matrix  equations 

If  -  Ye  and  M "  x  Y  M  (82) 

•9 

where  the  matrices  ME  and  MM  are  interface  matrices  for  TM  and  TE  waves,  respectively.  In 
accordance  with  the  notation  given  in  Section  3.3,  the  vectors  of  unknown  multipole 
coefficients  have  the  following  form: 


v  ^  ..  /mm  1  om  -am  J-jm  -am  -Ijm  — /im 

"ABe*  .v  di.r  d i.r  ^2.r  ^2,r  d2,r  d2,r  * 

\  pN,r  pN.r  f  V 


where  Xj  denotes  transpose  of  the  vector. 
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It  should  be  recalled  that  there  are  2  x  (2N  +  1)  elements  in  each  of  these  vectors.  The  form 
of  the  interface  matrices  is  shown,  schematically,  in  Fig.  4.  For  TM  waves,  the  rectangular 
submatrices  in  the  interface  matrix  have  the  following  form: 
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Figure  4.  Block  form  of  N  damper  interface  matrix  for  either 
TM  or  TE  waves. 
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Note  that  the  intervening  medium  is  not  required  to  be  free  space  but  only  nondissipative. 
The  submatrices  for  TE  waves  are  stated  below: 
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The  modular  nature  of  the  interface  matrices  ME  and  MM  is  dearly  demonstrated  in  Fig. 
4.  An  empty  cavity  is  represented  by  the  submatrices  M0,  Mex,  and  Mv  For  one  damper  the 
interface  matrices  are  built  up  from  M0,  MDp  Mex,  and  Ms.  A  cavity  containing  two 
dampers  is  represented  by  the  combination  MOJ  Md1?  MinI|,  MeK,  and  M,,.  A  three-damper 
cavity  is  represented  by  Mot  MDt,  Mimi,  MD2,  Min,2,  MD3,  Mex,  and  Ms,  etc.  Thus  a  damper 
can  be  added  to  or  subtracted  from  the  computational  scheme  merely  by  adding  or 
subtracting  a  damper  submatrix  MDj  and  a  submatrix  Mjnij,  corresponding  to  the 
intermediate  region  between  the  “j th”  damper  and  the  “j  +  \th"  damper. 

In  practice,  the  surface  conductivity  is  sufficiently  high  that  the  asymptotic  approxima¬ 
tions, 

and 

4  (jSsas)  3  el,lk«k"ff/2<f,+  ,)1i  =  (l  1  n/cajl 

can  be  used  to  simplify  the  corresponding  elements  in  Eqs.  (84e)  and  (85e).  No  such 
simplification  can  be  applied  to  the  remaining  matrix  elements,  since  a*  =  2.74  corresponds 
to  the  primary  TM  mode  of  any  spherical  cavity.  Nevertheless,  these  matrix  elements  can  be 
stably  computed  by  using  standard,  downward  recursion  algorithms  with  less  memory 
requirement  than  the  current  AEDC  Mie-scattering  codes  require. 

it  is  apparent  that  the  interface  matrices  are  sparse,  tetra-diagonal  matrices.  While  the 
innermost  two  diagonals  contain  exclusively  nonzero  elements,  the  outermost  two  contain 
alternating  zero  and  nonzero  elements.  This  structure  can  be  used  advantageously  to  devise 
efficient  inversion  algorithms  to  solve  Eq.  (81). 

There  are  three  significant  advantages  to  using  a  matrix  approach  to  solve  the  boundary 
value  problem  instead  of  using  the  direct  solutions  presented  in  the  next  section:  (1)  Because 
of  the  modular  nature  of  the  matrix,  the  number  of  dampers  considered  can  be  altered  by 
merely  adding  or  deleting  blocks  in  the  interface  matrix,  as  explained  earlier.  (2)  Less 
computational  difficulty  (and  less  program  complexity)  should  result  from  computing  the 
spherical  functions  used  in  Eqs.  (84)  and  (85)  (once,  for  each  frequency  of  interest)  and 
numerically  inverting  Eq.  (82)  to  obtain  the  desired  multipole  coefficients  (instead  of  using 
the  complicated  combinations  of  these  functions  appearing  in  the  next  section).  (3)  If  only 


(ps  a2N  +  l)\ 

<**> 
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the  frequency  response  of  the  source-cavity  system  (hence,  the  frequency  dependence  of  the 
antenna  impedance)  is  desired,  rather  than  the  electric  and  magnetic  field  distributions  in  the 
cavity,  then  only  the  standing  wave  coefficients  A™  and  Bjp  need  be  computed.  (Only  A" 
need  be  computed  for  an  axisymmetric  source.)  The  direct  solution  procedure,  in  effect, 
requires  that  all  other  multipole  coefficients  be  computed  in  the  process  of  determination  of 
A™  and  B£*.  Of  course  the  matrix  procedure  is  capable  of  simultaneously  determining  all 
other  multipole  coefficients,  but  they  need  not  be  computed,  unless  desired. 


3.3.2. 1  Multipole  Coefficients  for  a  Cavity  Containing  a  Single  Damper 


To  illustrate  the  direct  solution  procedure,  let  us  solve  the  single  damper  problem 
explicitly.  The  equations  that  relate  the  interior,  damped,  and  exterior  solutions  are 
analogous  to  Eq.  (75)  and  are  stated  below  for  TM  waves: 


The  simplest  way  to  solve  these  equations  is  to  work  backward  from  the  surface,  obtaining 
each  TM  wave  multipole  coefficient  in  terms  of  the  coefficient  in  the  immediately  exterior, 
adjacent  region.  This  very  cumbersome  procedure  yields  the  following  equations  [using  the 
asymptotic  forms  of  (Pgaa)  and  ffr  (/3sas)] : 
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The  coefficients  with  Roman  subscripts  in  Eq.  (87)  are  not  constants  but  are  the  auxiliary 
functions  stated  below: 


\ 
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(Note  that  the  Roman  subscripts  correspond  to  regions  in  Fig.  3.)  The  inconvenience  of 
using  these  direct  solutions,  in  contrast  to  the  apparent  simplicity  of  the  numerical  inversion 
procedure  advocated  in  the  previous  section,  is  annoyingly  obvious! 

The  TE  wave  multipole  coefficients  for  the  single-damper  configuration  are  obtained  in  a 
similar  fashion  and  are  stated  below: 
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The  associated  auxiliary  functions  are  the  following: 
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Since  most  treatments  of  damping  grids  assume  that  the  grids  are  infinitely  thin  surfaces 
across  which  the  tangential  magnetic  field  is  discontinuous  (with  this  discontinuity 
proportional  to  the  surface  current),  it  is  important  to  determine  when  the  finite  width  of  the 
dampers  must  be  considered.  Let  the  damper  width  be  D  =  rb  -  ra.  Now,  simplified  shifted 
argument  expansions  of  the  Ricatti-Hankel  functions  and  their  derivatives  can  be  stated 
(Ref.  5)  as 
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where  =  aj  -  or0  =  kD,  and  analogous  relations  can  be  derived  for  the  Ricatti-Bessel 
functions.  Equation  (91)  relates  the  damped  spherical  waves  at  the  inner  damper  radius  to 
those  at  the  outer  damper  radius.  After  some  manipulation  involving  the  Wronkskian,  W[f  * 
(0Dai)»  T?  (0d«i]  (see  Ref.  5),  we  find  that,  in  order  to  ignore  the  finite  width  of  the  damper 
in  the  computation  of  Cn  and  Dn>  the  following  inequalities  must  be  simultaneously 
satisfied: 
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For  the  lower-order  multipoles,  Eq.  (92)  simply  states  that  the  finite  width  of  the  damper  can 
be  ignored  if  it  is  much  less  than  Xd/2tt  where  XD  is  the  wavelength  of  the  TM  wave  in  the 
damper.  The  higher-order  (f  >  1)  multipoles,  however,  contribute  to  the  high  frequency 
response  of  the  system  because  Bessel  functions  peak  at  values  of  their  arguments  which  are 
comparable  to  their  order.  To  understand  the  significance  of  Eq.  (92),  we  consider  a  case 
where  the  damper  impedance  is  l/D,/2~times  that  of  the  intervening  medium  and  the  damper 
is  located  at  half  the  chamber  radius:  Then,  a,  =  0.7  as  where  ^/(as  fii)  =  0  corresponds 
to  the  TM  eigenmode  of  Ith  order  and  \th  degeneracy.  Consider  the  case  of  f  =  3,  since  even- 
order  multipoles  vanish  for  a  source  with  axisymmetry.  Then,  as  3>|  s  1.5  as  \y,  thus,  Eq. 
(92)  becomes 
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The  finite  width  of  the  damper  is  more  significant  when  the  chamber  resonance  is  of  a  high 
order  and  the  damper  is  near  the  center  of  the  chamber. 


If  desired,  the  preceding  Eqs.  (86)  through  (89)  can  be  simplified  by  retention  of  only 
first-order  terms  in  the  damper  width.  The  auxiliary  functions,  Cn  and  Dn,  reduce  to 
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with  similar  reductions  for  fn  and  gu.  Equations  (87d)  through  (87f)  are  then  evaluated  with 
a0  replaced  by  ai  and  using  Eq.  (88a)  for  Cn  and  D».  Obviously,  an  analogous  procedure 
can  be  applied  to  Eqs.  (89d)  through  (89f>- 

Simplifications  such  as  Eq.  (91)  can  also  be  applied  to  the  submatrices  in  Eqs.  (84)  and 
(85).  These  approximations  offer  the  obvious  advantage  of  reducing  the  number  of  spherical 
functions  to  be  computed.  In  addition,  they  isolate  the  small  parameter,  aD,  in  case  special 
precautions  must  be  taken  in  the  inversion  procedure  to  avoid  loss  of  significant  digits 
accuracy. 

3.3.2.2  Multipole  Coefficients  for  a  Cavity  Containing  Two  Dampers  - 

The  TM  wave  multipole  expansion  coefficients  (using  the  notation  convention  of  Section 
3.3)  for  the  two-damper  cavity  in  Fig.  5  are  given  below  (assuming,  for  simplicity,  that  the 
electrical  properties  of  both  dampers  are  the  same): 


Figure  5.  Geometry  for  electromagnetic  wave  propagation  in  an  internally 
excited  hollow  sphere  containing  two  concentric  damping  regions. 
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In  these  equations,  +7jn  are  the  multipole  coefficients  of 
outward  propagating  waves  at  the  inner  radius  of  the  sphere 
surface;  ej,",  inward  propagating  waves  at  the  inner  surface 
of  the  sphere;  Ef",  standing  waves  between  the  outer  damper 
and  the  sphere  inner  surface;  ± outward  and  inward 
propagating  waves,  respectively,  in  the  outer  damper;  +C1,1}, 
outward  and  inward  propagating  waves,  respectively,  in  the 
region  between  the  two  dampers;  outward  and  inward 

propagating  waves,  respectively,  in  the  inner  damper;  and 
Af",  standing  waves  between  the  inner  damper  and  the  center 
of  the  sphere  (exclusive  of  the  source). 
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There  are,  of  course,  analogous  equations  for  the  TE  wave  multipole  coefficients  and  their 
associated  auxiliary  functions,  but  these  will  not  be  stated. 

3.3.3  Coupling  of  the  Damped  Cavity  and  the  Biconical  Source 

As  shown  in  previous  sections,  the  multipole  coefficients  of  standing  waves  in  the  center 
of  a  damped  cavity  (cavity  multipole  coefficients)  are  related  to  the  source  multipole 
coefficients  through  equations  that  account  for  the  electrical  properties  of  all  other  regions 
of  the  sphere'.  Thus,  to  the  extent  that  linear  superposition  is  valid,  it  is  possible  to  calculate 
the  electric  and  magnetic  fields  in  the  vicinity  of  the  antenna  by  merely  adding  together  the 
contributions  of  the  source  multipoles  and  the  cavity  multipoles.  Strictly  speaking,  however, 
the  problem  is  not  linear,  since  the  presence  of  the  cavity  modifies  the  current  distribution  in 
the  source.  In  the  treatment  of  Sections  3.2. 1  and  3.2.2  this  interaction  has  been  neglected  in 
determining  the  source  current  distribution  and  the  associated  source  multipole  coefficients. 
Since  the  problem  is  virtually  intractable  if  this  interaction  is  included,  the  conditions  for 
which  the  linear  approximation  is  adequate  should  be  ascertained. 

Essentially  the  source  cannot  “see”  the  first  damper  until  the  leading  edge  of  the 
transmitted  pulse  reflects  from  the  damper  inner  surface  and  arrives,  again,  at  the  antenna. 
Thus,  linearity  is  an  adequate  assumption  for  all  times  t  <  (2ra/C).  Accordingly,  the  linear 
treatment  is  valid  for  all  circular  frequencies,  «  >  (C/2ra).  This  means  that  the  theory  is 
correct  for  reduced  damper  radii,  kra  >  1/2,  or  reduced  cavity  radii  (i.e.,  sphere  “size 
parameters”  in  optical  terminology),  kR  >  R/2ra).  Since  the  primary  TM  resonant  mode  of 
the  cavity  occurs  at  kR  =  2.74,  linear  coupling  is  valid  as  long  as  ra/R  >  (1/5.48)  =  0. 182, 
and  this  limit  has,  indeed,  not  been  violated  in  the  VKF  12-ft  sphere  experiments  in  which 
ra/R  =  0.8.  Obviously,  the  above  limit  must  not  be  disregarded  in  experimental  design, 
since  the  presence  of  any  conducting  surface  as  close  as  R/5.48  to  the  source  will  modify  the 
source  current  distribution  (at  frequencies  of  interest)  in  comparison  to  the  current 
distribution  calculated  when  it  is  assumed  that  the  source  radiates  into  an  unbounded 
medium. 

Assuming  the  validity  of  linear  superposition,  and  evaluating  Eq.  (17),  we  determine  the 
impedance  of  the  biconical  source  plus  damped  cavity  .to  be 

1)*  Pg(CO80e)  [a^'(l3>  +  Af  ^<k«] 

V  2C  +  1)54  Pt{c0sd^  [ae^1  <k*>+ 


ZA  =  * 
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where  Ar  =  A|?,  and  i(  should  be  recalled  that  only  odd  values  of  ( contribute  to  the  sums  in 
Eq.  (95),  on  account  of  the  cylindrical  symmetry  of  the  biconical  source.  The  frequency 
response  of  the  system  is,  then,  obtained  by  inserting  Eq.  (95)  into  Eq.  (19)  and  evaluating 
(numerically)  the  results  at  all  frequencies  greater  than  !comin  |  =  C/2ra. 

4.0  RECOMMENDATIONS  FOR  FUTURE  WORK 
AND  CONCLUDING  REMARKS 

This  report  has  presented  a  theoretical  formulation  to  describe  the  propagation  of  both 
TM  and  TE  waves  through  an  internally  excited,  damped  spherical  cavity.  The  general 
treatment  developed  herein  allows  consideration  of  an  arbitrary  central  source  current 
distribution  and  an  arbitrary  number  of  dampers  of  finite  width. 

Several  features  of  the  present  treatment  are  unique  and  worthy  of  further  comment: 

1.  The  present  theory  considers  dampers  with  distributed  electrical  properties  and 
finite  width,  and,  if  desired,  each  effective  damper  width  can  be  treated  as  an 
adjustable  parameter,  in  order  to  compensate  for  inhomogeneities  in  the  damping 
mesh.  Most  similar  investigations  have  assumed  that  the  dampers  can  be  treated  as 
infinitely  thin  surfaces  of  discontinuity.  While  the  finite  damper  width  is  not  very 
significant  in  analysis  of  the  12-ft  sphere  experiments,  some  future  analyses  of 
damping  configurations  may  require  inclusion  of  finite  width  effects. 

2.  The  present  theory  uses  a  matrix  formulation  to  obtain  the  multipole  expansion 
coefficients  throughout  the  regions  of  a  cavity  containing  N  dampers.  The 
computational  flexibility  afforded  by  the  modular  structure  of  this  matrix  permits 
easy  addition  or  deletion  of  dampers.  In  addition,  such  an  approach  is 
computationally  more  convenient  and  probably  more  stable  than  more  direct 
solution  procedures.  Although  this  matrix  approach  will  be  adopted  in  encoding 
the  theory,  explicit  solutions  for  the  single-damper  and  dual-damper  multipole 
expansion  coefficients  are  also  presented  in  this  report.  In  addition,  the  conditions 
under  which  the  finite  width  of  the  dampers  must  be  considered  are  investigated. 

3.  The  present  theory  uses  two  scalar  Hertz-Debye  potentials  to  describe  electric  and 
magnetic  fields  outside  the  source  region,  whereas  the  conventional  vector  and 
scalar  potentials  are  used  by  most  authors  There  are  two  advantages  associated  with 
the  use  of  these  potentials: 
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a.  The  propagating  fields  are  automatically  separated  into  transverse  magnetic 
(TM)  and  transverse  electric  (TE)  waves  in  any  geometry  in  which  Laplace’s 
equation  is  separable,  and 

b.  The  use  of  these  potentials  facilitates  the  solution  of  the  inhomogeneous 
wave  equations  which  represent  the  radiation  from  a  centralized,  arbitrary 
distribution  of  source  currents.  The  Hertz-Debye  potentials  relate  all  other 
electric  and  magnetic  field  components  to  the  radial  field  components,  whose 
solutions  are  obtained  by  the  use  of  Green’s  functions.  These  solutions  show 
that  only  TM  waves  are  radiated  by  cylindrically  symmetric  source  currents, 
and  TE  waves  are  radiated  either  if  an  azimuthal  current  flow  occurs  or  if 
azimuthal  gradients  in  polar  currents  occur. 

4.  The  present  theory  includes  a  closed-form  solution  for  the  multipole  expansion 
coefficients  which  represent  radiation  into  an  unbounded  medium  by  an 
uncapped  biconical  antenna.  In  addition,  the  two  coupled  integrodifferential 
equations— which  separately  describe  the  current  distribution  in  the  conical 
sections  and  in  the  caps  of  a  capped  biconical  antenna,  radiating  into  an 
unbounded  medium — are  derived  by  use  of  a  gauge  transformation.  These 
equations  can  be  solved  by  the  method  of  moments  or  by  a  similar  finite-element 
numerical  technique.  The  multipole  expansion  coefficients  for  a  capped  bicone 
can  then  be  obtained  by  substituting  the  bicone  current  expansion  into  the  TM 
wave  Green’s  function  solution. 

5.  The  frequency  response  (and,  by  inverse  Fourier  transformation,  the  time 
dependent  pulse  waveform)  of  the  internally  excited,  damped  cavity  is  obtained 
from  the  characteristic  impedance  of  the  bicone  and  the  input  impedance  of  the 
antenna,  radiating  into  the  damped  cavity.  This  impedance,  in  turn,  is  calculated 
from  the  multipole  expansion  coefficients  of  the  source  and  the  first  standing 
waves  in  the  region  between  the  source  and  the  damper.  The  antenna  impedance 
and  frequency  response  of  the  antenna,  radiating  into  an  unbounded  medium, 
can  be  obtained  simply  by  setting  the  multipole  coefficients  of  the  standing 
waves  equal  to  zero.  This  procedure  provides  a  check  on  the  accuracy  of  the 
source  current  solution. 

Future  work  in  this  area  should  follow  the  directions  indicated  below: 

1.  Encode  the  theory,  including  the  bicone  source  current  determination,  with 
options  to  calculate  the  (polar)  angular  and  radial  electric  and  magnetic  field 
distributions  throughout  the  cavity  at  specified  frequencies  and  to  calculate  only 
the  frequency  response  of  the  cavity  at  the  bicone  feed  terminals  (hence,  by 
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inverse  Fourier  transformation,  the  wave  form  of  the  sensor  current  pulse).  The 
primary  result  of  this  work  would  be  a  comparison  between  theory  and 
experiment.  The  tasks  that  follow  cannot  be  implemented  until  this  has  been 
accomplished. 

2.  Adapt  the  theory  and  its  coding  to  cylindrical  coordinates. 

3.  Account  for  the  effect  of  openings  to  accommodate  instrumentation  and  X-ray 
beam  paths  in  the  damping  mesh.  In  an  infinitely  thin  damper,  such  openings  are 
denoted  “Bethe  appertures”  when  they  are  smaller  than  a  wavelength,  and  they 
act  as  dipole  sources  located  at  the  positions  of  the  openings. 

4.  Generalize  the  theory  to  account  for  an  off-centered,  but  still  axisymmetric 
source. 
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NOMENCLATURE 

A™  Standing  TM  wave  multipole  coefficients  for  central  cavity  region 

(exclusive  of  source) 

a™  Source  TM  wave  multipole  coefficients 

B  Magnetic  induction 

B™  Standing  TE  wave  multipole  coefficients  for  central  cavity  region 

(exclusive  of  source) 

b™  Source  TE  wave  multipole  coefficients 

Cn .  v  TM  wave  auxiliary  functions  defined  in  Eqs.  (88)  and  (94) 

c  Speed  of  light  in  vacuum 

±  c™  Outwardly  and  inwardly  propagating  TM  wave  multipole  coefficients  for 

region  between  “j th”  and  “j  +  Ith"  dampers 

D  Damper  thickness 

D  Displacement 

Du  .  v  Derivatives  of  TM  wave  auxiliary  functions  defined  in  Eqs.  (88)  and  (94) 

±  d™  Outwardly  and  inwardly  propagating  TE  wave  multipole  coefficients  for 

j  region  between  “j  th  and  “j  +  Uh"  dampers 

E  Electric  field 

EJP  Standing  TM  wave  multipole  coefficients  for  region  between  outermost 

damper  and  sphere  surface 

\ 

±  e™  Outwardly  and  inwardly  propagating  TM  wave  multipole  coefficients  for 

region  between  outermost  damper  and  sphere  surface 

fn  m  TE  wave  auxiliary  functions  defined  in  Eq.  (90) 
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Green’s  functions  used  in  solution  of  inhomogeneous  Helmholtz 
equations 

Derivatives  of  TE  wave  auxiliary  functions  defined  in  Eq.  (90) 

Magnetic  field 

Axial  length  of  single  cone  of  a  biconical  antenna  (see  Fig.  2) 

l 

Electric  current 

Complex  constant  appearing  in  Eqs.  (40)  through  (42) 

Electric  current  density 
Integral  defined  in  Eq.  (45) 

Source  current  density 

Kernel  function  defined  in  Eq.  (71) 

Wave  propagation  constant 

Wave  propagation  constant  for  a  nondissipative  medium 
Radial  length  of  single  cone  of  a  biconical  antenna  (see  Fig.  2) 

Angular  momentum  operator 
Kernel  function  defined  in  Eq.  (71) 

Interface  matrices  defined  in  Eqs.  (84)  and  (85) 

Standing  TE  wave  multipole  coefficients  for  region  between  outermost 
damper  and  sphere  surface 

Inwardly  propagating  TE  wave  multipole  coefficients  for  region  between 
outermost  damper  and  sphere  surface 

Kernel  function  defined  in  Eq.  (71) 
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Pe,m 
Pc  (cos  0) 
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Ru 
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TE 
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Number  of  dampers 

Hertz  vectors  defined  by  Eq.  (19) 

Legendre  polynomial  of  order  ( 

Sphere  radius 

Distances  defined  in  Eqs.  (56)  and  (57) 

Position  vector  to  field  observation  point 

Position  vector  to  source  point 

Scalar  function  defined  in  conjunction  with  Eq.  (25) 

Arbitrary  scalar  function 
Transverse  electric 
Transverse  magnetic 
Time 

Thickness  of  current-carrying  region  in  sides  of  antenna 
Thickness  of  current-carrying  region  in  caps  of  antenna 
Transformed  field  defined  in  Eq.  (49) 

Potential  difference  between  the  two  halves  of  a  biconical  antenna 
Arbitrary  vector  function 

Column  vectors  of  unknown  multipole  coefficients  [see  Eq.  (83)] 
Column  vectors  associated  with  Eq.  (83) 

Spherical  harmonic  functions 
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Wave  impedance 

i 

Characteristic  impedance  of  a  beconical  antenna 
Impedance  of  a  nondissipative  medium 
Impedance  of  free  space 
Coordinate  along  axis  of  symmetry 
Dimensionless  radius  of  the  “j  th"  interface 

Outwardly  and  inwardly  propagating  TM  wave  multipole  coefficients  at 
the  interface 

Ratio  of  propagation  constants  in  dissipative  and  nondissipative  media 

Outwardly  and  inwardly  propagating  TE  wave  multipole  coefficients  at 
the  “jth"  interface 

Outwardly  propagating  TM  wave  multipole  coefficients  at  the  sphere 
surface 

Radial  location  of  antenna  tip  {see  Fig.  2) 

Penetration  depth 

Outwardly  propagating  TE  wave  multipole  coefficients  at  the  sphere 
surface 

Real  dielectric  permittivity 
Complex  dielectric  permittivity 
Electrical  conductivity 
Magnetic  permeability 

Hertz-Debye  potentials  for  TM  and  TE  waves 
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Complementary  antenna  function  defined  in  Eq.  (48) 

Modified  Green’s  functions  defined  in  Eq.  (64) 

Antenna  cylindrical  coordinate 
Charge  density 
Radius  of  antenna  cap 
Gauge  variable  defined  in  Eq.  (49) 

Ricatti-Bessel  functions  of  order  ( 

Ricatti-Hankel  functions  or  order  f 
Wavelength  , 

Circular  frequency 

SUPERSCRIPTS 

E  Pertains  to  electric  Hertz-Debye  potential  (source  of  TM  waves) 

M  Pertains  to  magnetic  Hertz-Debye  potential  (source  of  TE  waves) 

m  Azimuthal  index  for  spherical  expansion 

T  Denotes  matrix  transpose 

Denotes  source  coordinate  (except  in  Sec.  3.3,  where  it  denotes 
differentiation) 

~  Denotes  frequency  domain 

+  Denotes  outwardly  propagating  wave 


Xr(k»kpid) 

xtftfr'.W 
xtftfr',*') 
x I* (r ,  e.*') 

Xfc  G. 
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Qc 
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-  Denotes  inwardly  propagating  wave 

Denotes  unit  vector 
-**  Denotes  vector 

SUBSCRIPTS 

A  Pertains  to  antenna 

be  Pertains  to  conical  surfaces  of  biconical  antenna 

cp  Pertains  to  cap  surfaces  of  biconical  antenna 

D  Pertains  to  damper 

E  Pertains  to  electric  Hertz-Debye  potential  (source  of  TM  waves) 

II,  III,  ...  Pertain  to  cavity  regions  shown  in  Figs.  3  and  5 
f  Polar  index  for  spherical  expansion 

M  Pertains  to  magnetic  Hertz-Debye  potential  (source  of  TE  waves) 

s  Denotes  quantity  evaluated  at  surface  of  cavity 

_  Denotes  matrix  or  column  vector 
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